Bicritical Behavior of Period Doublings in Unidirectionally-Coupled Maps 
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We study the scaling beliavior of period doublings in two unidirectionally-coupled one- 
dimensional maps near a bicritical point where two critical lines of period-doubling transition to 
chaos in both subsystems meet. Note that the bicritical point corresponds to a border of chaos 
in both subsystems. For this bicritical case, the second response subsystem exhibits a new type 
of non-Feigenbaum critical behavior, while the first drive subsystem is in the Feigenbaum critical 
state. Using two different methods, we make the renormalization group analysis of the bicritical 
behavior and find the corresponding fixed point of the renormalization transformation with two 
relevant eigenvalues. The scaling factors obtained by the renormalization group analysis agree well 
with those obtained by a direct numerical method. 

PACS numbers: 05.45. +b, 03.20.-|-i, 05.70.Jk 



I. INTRODUCTION 

Period-doubling transition to chaos has been ex- 
tensively studied in a one-parameter family of one- 
dimensional (ID) unimodal maps, 



Xt+l 



1 - Axf 



(1) 



where Xt is a state variable at a discrete time t. As the 
control parameter A is increased, the ID map undergoes 
an infinite sequence of period-doubling bifurcations accu- 
mulating at a critical point Ac, beyond which chaos sets 
in. Using a renormalization group (RG) method, Feigen- 
baum Q has discovered universal scaling behavior near 
the critical point Ac- 
Here we are interested in the period doublings in a sys- 
tem consisting of two ID maps with a one-way coupling. 



xt+i = 1 - Axf , yt+i = 1- By, 



Cxi , 



(2) 



where x and y are state variables of the first and sec- 
ond subsystems, A and B are control parameters of the 
subsystems, and C is a coupling parameter. Note that 
the first (drive) subsystem acts on the second (response) 
subsystem, while the second subsystem does not influ- 
ence the first subsystem. This kind of unidirectionally- 
coupled ID maps are used as a model for open flow |H . A 
new kind of non-Feigenbaum scaling behavior was found 
in a numerical and empirical way near a bicritical point 
{Ac, Be) where two critical lines of period-doubling tran- 
sition to chaos in both subsystems meet For this 
bicritical RG analysis was also developed and the 

corresponding fixed point, governing the bicritical behav- 
ior, was numerically obtained by directly solving the RG 
fixed-point equation using a polynomial approximation 



Q. In this paper, using two different methods, we also 
make the RG analysis of the bicriticality, the results of 
which agree well with those of previous works. 

This paper is organized as follows. In Sec. || we study 
the scaling behavior near a bicritical point {Ac, Be), cor- 
responding to a border of chaos in both subsystems, by 
directly following a period-doubling sequence converging 
to the point {Ac, Be) for a fixed value of C. For this 
bicritical case, a new type of non-Feigenbaum critical be- 
havior appears in the second subsystem, while the first 
subsystem is in the Feigenbaum critical state. Employing 
two different methods, we make the RG analysis of the bi- 
critical behavior in Sec. m. To solve the RG fixed-point 
equation, we first use an approximate truncation method 

, corresponding to the lowest-order polynomial approx- 
imation. Thus we analytically obtain the fixed point, 
associated with the bicritical behavior, and its relevant 
eigenvalues. Compared with the previous numerical re- 
sults [Q , these analytic results are not bad as the lowest- 
order approximation. To improve accuracy, we also em- 
ploy the "eigenvalue-matching" RG method j|] , equating 
the stability multipliers of the orbit of level n (period 2") 
to those of the orbit of the next level n -|- 1. Thus we nu- 
merically obtain the bicritical point, the parameter and 
orbital scaling factors, and the critical stability multipli- 
ers. We note that the accuracy is improved remarkably 
with increasing the level n. Finally, a summary is given 
in Sec. fv. 



II. 



SCALING BEHAVIOR NEAR THE 
BICRITICAL POINT 



In this section we fix the value of the coupling pa- 
rameter by setting C = 0.45 and directly follow a 
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period-doubling sequence converging to the bicritical 
point {Ac-, Be), which corresponds to a border of chaos 
in both subsystems. For this bicritical case, the second 
subsystem exhibits a new type of non-Feigenbaum critical 
behavior, while the first subsystem is in the Feigenbaum 
critical state. 

The unidirectionally-coupled ID maps (Q) has many 
attractors for fixed values of the parameters . For the 
case C = 0, it breaks up into the two uncoupled ID maps. 
If they both have stable orbits of period 2*^ , then the com- 
posite system has 2*^ different stable states distinguished 
by the phase shift between the subsystems. This mul- 
tistability is preserved when the coupling is introduced, 
at least while its value is small enough. Here we study 
only the attractors whose basins include the origin (0, 0). 
Such attractors become in-phase when A = B and C = 0. 

Stability of an orbit with period q is determined by its 
stability multipliers. 



Ai = fl -2Axt 



n -^Byt 



(3) 



Here Ai and A2 determine the stability of the first and 
second subsystems, respectively. An orbit becomes stable 
when the moduli of both multipliers are less than unity, 
i.e., -1 < Ai < 1 for i = 1,2. 



first subsystem exhibits a sequence of period-doubling bi- 
furcations at the vertical straight lines, where Ai = — 1. 
For small values of the parameter B, the period of os- 
cillation in the second subsystem is the same as that in 
the first subsystem, as in the case of forced oscillation. 
As B is increased for a fixed value of A, a sequence of 
period-doubling bifurcations occurs in the second subsys- 
tem when crossing the non- vertical lines where A2 = —1. 
The numbers inside the different regions denote the pe- 
riod of the oscillation in the second subsystem. 

We consider a pair of the parameters (^„,i?„), at 
which the periodic orbit of level n (period 2") has the 
stability multipliers Ai_„ — A2,n = —1- Hence, the point 
{An, Bn) corresponds to a threshold of instability in both 
subsystems. Some of such points are denoted by the open 
circles in Fig. 1. Then such a sequence of {An,Bn) con- 
verges to the bicritical point {Ac, Be), corresponding to 
a border of chaos in both systems, with increasing the 
level n. The bicritical point is denoted by the solid circle 
in Fig. 1. To locate the bicritical point with a satisfac- 
tory precision, we numerically follow the orbits of period 
g = 2" up to level n = 21 in a quadruple precision, and 
obtain the sequences of both the parameters (^„,_B„) 
and the orbit points {xn, Vn) approaching the origin. We 
first note that the sequences of An and Xn in the first sub- 
system are the same as those in the ID maps Hence, 
only the sequences of Bn and yn in the second subsystem 
are given in Table |. 
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TABLE I. 


Sequences of the parameter 


and the orbit point. 


{Bn} and {yn}, in the second subsystem. 






n 




yn 




10 


1.090 088 955 364 


5.019189 X 10" 


-■J 


11 


1.090 092 109 910 


-3.333 775 xlO 


-3 


12 


1.090 093 416 851 


2.214 467 X 10" 


-3 


13 


1.090 093 959 979 


-1.471024 X 10 


-3 


14 


1.090 094186 392 


9.771 970 xlO" 


-4 


15 


1.090 094 280 906 


-6.491561 XlO 


-4 


16 


1.090 094 320 376 


4.312 391 XlO" 


-4 


17 


1.090 094 336 865 


-2.864 762 xlO 


-4 


18 


1.090 094 343 755 


1.903 092 XlO" 


-4 


19 


1.090 094 346 634 


-1.264 245 XlO 


-4 


20 


1.090 094 347 837 


8.398 518 XlO" 


-5 


21 


1.090 094 348 340 


-5.579 230 XlO 


-5 



1.36 



1.38 



1.40 



FIG. 1. Stability diagram of the periodic orbits born via 
period-doubling bifurcations for C = 0.45. The numbers in 
the different regions represent the period of motion in the sec- 
ond subsystem. The open circle also denotes the point, cor- 
responding to a threshold of instability in both subsystems, 
where Ai = —1 and A2 — —1. Such open circles accumulate 
to the bicritical point, denoted by the solid circle, which cor- 
responds to a border of chaos in both subsystems. For other 
details, see the text. 

Figure |l| shows the stability diagram of periodic or- 
bits for C = 0.45. As the parameter A is increased, the 



We now study the asymptotic scaling behavior of the 
period-doubling sequences in both subsystems near the 
bicritical point. The scaling behavior in the first sub- 
system is obviously the same as that in the ID maps 
[§. That is, the sequences {An} and {a;„} accumulate to 
their limit values, A ^ Ac {^ 1.401 155 189 092 •• •) and 
X = 0, geometrically as follows: 



An Ac ^ , 



for large n. 



(4) 



The scaling factors Si and ai are just the Feigenbaum 
constants S {= 4.669 • • •) and a {= -2.502 • • •) for the ID 
maps, respectively. However, the second subsystem ex- 



2 



hibits a non-Feigenbaum critical behavior, unlike the case 
of the first subsystem. The two sequences {Bn} and {j/n} 
also converge geometrically to their limit values B — Be 
(= f.090 094 348 70f) and y = 0, respectively, where the 
value of Be is obtained using the superconverging method 
1^ . To obtain the convergence rates of the two sequences, 
we define the scaling factors of level n: 

r _ Bn-l — Bn _ Vn-l — Un ,r-\ 
02,n = ^ 5 , a2,n = • (5) 



Bn — Bn+1 



Un - Vn+l 



These two sequences {(52, «} and {a2.n} are listed in Table 
H, and they converge to their limit values. 



62 ^ 2.3928, a2 



-1.5053, 



(6) 



the region in the small box (containing the origin) by 
the scaling factor ai for the x axis and a2 for the y 
axis, and then we get the picture in Fig. |2|(c). Note that 
the picture in Fig. ^(c) reproduces the previous one in 
Fig. |^(a) approximately. Repeating the above procedure 
once more, wc obtain the two pictures in Figs. |^(d) and 
||(e). That is. Fig. ||(d) shows the hyperchaotic attractor 
with cr ~ 0.061 and era ^ 0.01 for A = Ae + AA/cJ^ and 
B = _Bc + A_B/(S|. Magnifying the region in the small box 
with the scaling factors a\ for the a;-axis and for the 
y-axis, we also obtain the picture in Fig. |^(e), which re- 
produces the previous one in Fig. |^(c) with an increased 
accuracy. 



respectively. Note that these scaling factors are com- 
pletely different from those in the first subsystem (i.e., 
the Feigenbaum constants for the ID maps). 

TABLE II. Sequences of the parameter and orbital scaling 
factors, {52, n} and {Q2,n}, in the second subsystem. 



02,; 



10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 



2.429 8 
2.413 7 
2.406 3 
2.398 8 
2.395 6 
2.394 6 
2.393 7 
2.393 1 
2.393 
2.392 9 
2.392 8 



-1.505 7331 
-1.505 515 4 
-1.505 4281 
-1.505 375 3 
-1.505 344 
-1.505 3316 
-1.505 325 6 
-1.505 3215 
-1.505 319 8 
-1.505 3191 
-1.505 318 6 



For evidence of scaling, we compare the chaotic attrac- 
tors, shown in Fig. ||, for the three values of {A^ B) near 
the bicritical point {Ac, Be). All these attractors are the 
hyperchaotic ones with two positive Lyapunov exponents 



ni ^ ni 

(71= lim — In |2^a;t|, (T2 = lim — y^ln|2_Bj/t 



(7) 



Here the first and second Lyapunov exponents a\ and 
02 denote the average exponential divergence rates of 
nearby orbits in the the first and second subsystems, re- 
spectively. Figure ||(a) shows the hyperchaotic attractor 
with cTi ~ 0.242 and 02 ^ 0.04 for A = Ac + AA and 
B = Bc + AS, where A A = AB = 0.1. This attrac- 
tor consists of two pieces. To see scaling, we first rescale 
AA and Ai? with the parameter scaling factors and 
62, respectively. The attractor for the rescaled parame- 
ter values of A = Ac + A A/ Si and B = Bc + AB/S2 is 
shown in Fig. |^(b). ft is also the hyperchaotic attrac- 
tor with (Ti ~ 0.121 and (T2 — 0.02. We next magnify 
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FIG. 2. Hyperchaotic attractors for the three val- 
ues of {A,B) near the bicritical point (Ac-Be); in (a) 
{A,B) = {Ac + AA,Bc + AB) {AA = AB = 0.1), in (b) 
and (c) {A,B) = {Ac + AA/Si, Be + AB /52), and in (d) and 
(e) {A, B) = {Ac + AA/Sf, Be + AB/Sl). The picture in (c) 
is obtained by magnifying the region in the small box in (b) 
with the scaling factors Qi for the a;-axis and 02 for the y-axis. 
Similarly, we also obtain the picutre (e) by magnifying the re- 
gion inside the small box in (d) with the scaling factors ol\ for 
the a;-axis and a\ for the j/-axis. Comparing the pictures in 
(a), (c), and (e), one can see that each successive magnified 
picture reproduces the previous one with an accuracy with 
the depth of resolution. 

So far we have seen the scaling near the bicritical point, 
and now turn to a discussion of the behavior exactly at 
the bicritical point {Ac^Bc). There exist an infinity of 
unstable periodic orbits with period 2" at the bicritical 



3 



point. The orbit points a;„ and y„, approaching the zero 
in the first and second subsystems, vary asymptotically 
in proportion to a^" and aj"", respectively. The stabil- 
ity multipliers Ai_„ and \2,n of the orbits with period 2" 
also converge to the critical stability multipliers and 
A2, respectively. Here A^ (= —1.601 191 • • •) in the first 
subsystem is just the critical stability multiplier for the 
case of the ID maps [Q. However, as listed in Table [H, 
the second subsystem has the different critical stability 
multiplier. 



A^ = -1.178 85 



(8) 



Consequently, the periodic orbits at the bicritical point 
have the same stability multipliers \\ and A2 for suffi- 
ciently large n. 

TABLE III. Sequences of the second stability multipliers, 
{A2,n} of the orbits with period 2" at the bicritical point. 

n Ao.n 



10 
11 
12 
13 
14 
15 
16 
17 
18 
19 



-1.178 829 
-1.178 842 
-1.178 839 
-1.178 850 
-1.178 855 
-1.178 854 
-1.178 854 
-1.178 855 
-1.178 855 
-1.178 854 



III. RENORMALIZATION GROUP ANALYSIS 
OF THE BICRITICAL BEHAVIOR 

Employing two different methods, we make the RG 
analysis of the bicritical behavior. We first use the trun- 
cation method, and analytically obtain the correspond- 
ing fixed point and its relevant eigenvalues. These ana- 
lytic results are not bad as the lowest-order approxima- 
tion. To improve the accuracy, we also use the numeri- 
cal eigenvalue-matching method, and obtain the bicritical 
point, the parameter and orbital scaling factors, and the 
critical stability multipliers. Note that the accuracy in 
the numerical RG results is improved remarkably with 
increasing the level n. 



A. Truncation Method 

In this subsection, employing the truncation method 
we analytically make the RG analysis of the bicriti- 
cal behavior in the unidirectionally-coupled map T of the 
form, 



where Xt and yt are the state variables at a discrete time 
t in the first and second subsystems, respectively. Trun- 
cating the map (0) at its quadratic terms, we have 



xt+i ^ - + bx 
b 



yt+i ^ -+dy. 



(10) 



which is a five-parameter family of unidirecionally- 
coupled maps. P represents the five parameters, i.e., 
P = (a, 5, c, c?, e). The construction of Eq. (^ corre- 
sponds to a truncation of the infinite dimensional space 
of unidirectionally-coupled maps to a five-dimensional 
space. The parameters a, 6, c, d, and e can be regarded as 
the coordinates of the truncated space. We also note that 
this truncation method corresponds to the lowest-order 
polynomial approximation. 

We look for fixed points of the renormalization op- 
erator TZ in the truncated five-dimensional space of 
unidirectionally-coupled maps. 



7^(T) = AT^A~\ 
Here the rescaling operator A is given by 



A 



ai 
a2 



(11) 



(12) 



where ai and a2 are the rescaling factors in the first and 
second subsystems, respectively. 

The operation TZ in the truncated space can be rep- 
resented by a transformation of parameters, i.e., a map 
from P = (a, 5, c, d, e) to P' = (a', 6', c', d' , e'), 



2a^{l + a), 

2 , 

— ao, 
ai 

2c(c + + 
— cd, 

0-2 

4 , 

-^ce[a -|- c). 



(13a) 
(13b) 

(13c) 
(13d) 
(13e) 



The fixed point P* = {a* ,b* ,c* , d* , e*) of this map can 
be determined by solving P' = P. The parameters b and 
d set only the scales in the x and y, respectively, and thus 
they are arbitrary. We now fix the scales in x an d y b y 
setting b — d = 1. Then, we have, from Eqs.(13a)-(13e), 
five equations for the five unknowns ai, a*, a2, c*, and e*. 
We thus find one solution, associated with the bicritical 
behavior, as will be seen below. The map (|l^ with a so- 
lution P* (Tp« ) is the fixed map of the renormalization 
transformation TZ; for brevit y Tp * will be denoted as T*. 

We first note that Eqs. ( |l3a| )-( p^ ) are only for the 
unknowns ai and a*. We find one solution for ai and 
a* , associated with the period-doubling bifurcation in the 
first subsystem. 



T : xt+i = f{xt), yt+i ^ gixt,yt), 



(9) 



ai 



= -1 - \/3 = -2.732 • • • , a* = Y 



(14) 
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Subs tituting the values for ai and a* into Eqs. (13c)- 
(13e), we obtain one solution for Q!2, c*, and e*, associ 
ated with the bicriticality, 



a2 



i(l + \/3-\/5-\/l5) = -1.688- ••, 



(15a) 



c* = e* = 1 + ^(VTS- SVS) = 0.338- • •. (15b) 

Compared with the values, ai ~ —2.502 • • • and a2 = 
— 1.505 • • •, obtained by a direct numerical method, the 
analyitc results for ai and a2, given in Eqs. ( p^ and 
(15a), are not bad as the lowest-order approximation. 



Consider an infinitesimal perturbation e SF to a fixed 
point P* of the transformation of parameters (13a)-(13e). 
Linearizing the transformation at P*, we obtain the 
equation for the evolution of SP, 



5P' = JSP, 



(16) 



where J is the Jacobian matrix of the transformation at 



The 5x5 Jacobian matrix J has a semi-block form, 
because we are considering the unidirectionally-coupled 
case. Therefore, one can easily obtain its eigenvalues. 
The first two eigenvalues, associated with the first sub- 
system, are those of the following 2x2 matrix, 



Ml = 



d{a',b') 



d{a,b) 



3-ai 
2/ai 1 



(17) 



Hence the two eigenvalues of Mi , 5i and 5'i , are given by 

(5i =4 + V3 = 5.732---, 5'i^l. (18) 

Here the relevant eigenvalue 5i is associated with the 
scaling of the control parameter in the first subsystem, 
while the marginal eigenvalue 5'i is associated with the 
scale change in x. When compared with the numerical 
value, ^1 (= 4.669 • • •), the analytic result for (5i, given in 
Eq. (Eq), is not bad as the lowest-order approximation. 

The remaining three eigenvalues, associated with the 
second subsystem, are those of the following 3x3 matrix. 



M2 = 



^(c^d^eo 

d(c, d, e) 



4c* + 6d 



,*2 



2a 



*2, 



2a 

2/a2 2c*/a2 
4:e*{a* + 2c*)/al 4c*(a 



c*)/al 



(19) 



The three eigenvalues of M2, 62, S'2, and 62, are given 

by 

^2 = (m -I- V^) /2 = 3.0246 • • • , (20a) 
62 = {u - Vv)/2 = 0.1379- ■■, 62=1, (20b) 

where 

u= (17 + 7\/3-5\/5-3%/T5)/2, (21a) 
t; = 104 + 53%/3-44V5-23VT5. (21b) 

The first eigenvalue 82 is a relevant eigenvalue, associated 
with the scaling of the control parameter in the second 
subsystem, the second eigenvalue 62 is an irrelavant one, 
and the third eigenvalue S2 is a marginal eigenvalue, as- 
sociated with the scale change in y. We also compare 
the analytic result for 62, given in Eq. ( |20a| ), with the 
numer ical value (^2 - 2.3928) in Eq. (|), and find that 
the analytic one is not bad as the lowest-order approxi- 
mation. 

As shown in Sec. stability multipliers of an orbit 
with period 2" at the bicritical point converge to the 
critical stability multipliers, (= —1.601 • • •) and 
(= —1.178 • • •) as n — > 00. We now obtain these critical 
stability analytically. The invariance of the fixed map 
T* under the renormalization transformation TZ implies 
that, if T* has a periodic point (x,y) with period 2", 
then A~^(a;, y) is a periodic point of T* with period 2"+^. 
Since rescaling does not affect the stability multipliers, all 



the orbits with period 2" {n = 0, 1, 2, . . .) have the same 
stability multipliers, which are just the critical stability 
multipliers A^ and A2. That is, the critical stability mul- 
tipliers have the values of the stability multipliers of the 
fixed point {x*,y*) of the fixed map T*, 

At = 2a;* = -1.5424- • •, A^ = 2?/* = -0.8899 ••• , (22) 

where 

X* = (1 - \/3 + 2V3)/2, y* - (1 - V^)/2, (23a) 

w^5 + 3\/3 - 2V5 - Vl5 

+ \/3 + 2V3 (2 - 3\/3 + VT5). (23b) 

We also note that the analytic values for A| and A2 are 
not bad, when compared with their numerical values. 

B. Eigenvalue-Matching Method 

In this subsection, we employ the eigenvalue-matching 
method Q and numerically make the RG analysis of the 
bicritical behavior in the unidirectionally-coupled map T 
of Eq. (^). As the level n increases, the accuracy in the 
numerical RG results are remarkably improved. 

The basic idea is to associate a value {A',B') for 

each value (A, B) such that T^2^^') locally resembles 
T^2^B), where T^") is the 2"th-iterated map of T (i.e.. 
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rp(n) =7^2") ^ simple way to implement this idea is to 
linearize the maps in the neighborhood of their respective 
fixed points and equate the corresponding eigenvalues. 

Let {zt} and {z^} be two successive cycles of period 2" 
and 2"+^, respectively, i.e., 



^* = ^(lB)(^t)' 4^T^2'.B')i^iy, ixt,yt)- (24) 

Here Xt depends only on A, but yt is dependent on both 
A and B, i.e., Xt = Xt{A) and yt — yt{A, B). Then their 
linearized maps at zt and z[ are given by 



DTliB)-X{mA,B){^t). 



t=l 

2"+i 



DT, 



(n+l) 



{A\B' 



n DT(A',B'M). 



(25a) 
(25b) 



(Here DT is the linearized map of T.) Let 
their eigenvalues, called the stability multipliers, be 
(Ai,„(A),A2,„(A,B)) and (Ai,„+i(A'), A2,„+i(A', S'))- 
The recurrence relations for the old and new parame- 
ters are then given by equating the stability multipliers 
of level n, Ai^„(A) and A2,ri(^, -B), to those of the next 
level n+l, Ai^„+i(A') and A2,„+i(A', S'), i.e.. 



Ai,„(A) = Ai,„+i(A'), 

A2,„(A^) = A2,„+l(A',B'). 



(26a) 
(26b) 



The fixed point (A*, B*) of the renormalization trans- 
formation (E6h, 



Ai,„(A*) = Ai,„+i(A*), 

A2,„(A*,B*) = A2,„+l(A*,B* 



(27a) 
(27b) 



gives the bicritical point {Ac, Be). By linearizing the 
renormalization transformation ( |2q ) at the fixed point 
(A*, B*), we have 







9A 




* 


dB' 


dB 




dB 


dA' 




dB' 


( AA 





' i AS 





A^' 
AS' 



(28) 
(29) 



where AA = A - A*, AS = B - B\ AA' = A' - A*, 
AS' = B' - B*, and 

(30a) 

r„ = I .r * I , (30b) 



Tri+i; 






dAi,„ 









dA 




.) 


dX2.n. 




aA2,„ 




dA 




dB 




c!Ai,„+i 







dA' 






9A2,„+i 


9A2 


n + l 


9A' 







(30c) 



Here F^^ is the inverse of r„ and the asterisk denotes 
the fixed point {A*,B*). After some algebra, we obtain 



the analytic formulas for the eigenvalues 5i.„ and S2.n of 
the matrix A„, 



S2.71 = 



dAi_„ + i 




dA' 




d\i,„ 






dA 





d\2 



dB' 



d\2 



(31a) 



(31b) 



As n — > cx) (Si.„ and S2^n approach Si and S2, which are 
just the parameter scaling factors in the first and second 
subsystems, respectively. Note also that as in the ID 
case, the local rescaling factors of the state variables are 
simply given by 



012.'. 



dx 
dx' 

dy 

dy' 



Sl,n 
tl,n 
^2,n 
^2,n 



(32a) 
(32b) 



where 



dx' 




dA' 








dx 




dA 


* 



dy' 

dB' 



(33) 



Here q;i.„ and a2,n also converge to the orbital scaling 
factors, cki and 012, in the first and second subsystems, 
respectively. 

We numerically follow the orbits with period 2" in the 
unidirectionally-coupled maps T of Eq. (|^) and make the 
RG analysis of the bicritical behavior. Some results for 
the intermediate level n are shown in Fig. ^. Figure |^(a) 
shows the plots of the first stability multiplier Ai^„(^) 
versus A for the cases n = 6, 7. We note that the in- 
tersection point, denoted by the solid circle, of the two 
curves Ai^g and Ai,7 gives the point (^g. A* g) of level 6, 
where Aq and A* g are the critical point and the critical 
stability multiplier, respectively, in the first subsystem. 
As the level n increases, A* and A* „ approach their limit 
values A* and A*, respectively. Note also that the ra- 
tio of the slopes of the curves, Ai_6(^) and Ai_7(^), for 
yl = ^g gives the parameter scaling factor <5i,6 of level 
6 in the first subsystem. Similarly, Fig. ||(b) shows the 
plots of the second stability multiplier A2,r! (v4g, B) versus 
B for the cases n — Q,l . The intersection point, denoted 
also by the solid circle, of the two curves A2,6(^6' ^) ^^d 
X2j{Al,B) gives the point (i3g,A2g) of level 6, where 
Bg and g are the critical point and the critical stabil- 
ity multiplier, respectively, in the second subsystem. As 
the level n increases, i?* and A2 „ also converge to their 
limit values, B* and A2, respectively. As in the first sub- 
system, the ratio of the slopes of the curves, A2,6(^6' ^) 
and A2,7(^g,i?), for B = B^ gives the parameter scaling 
factor 62,6 of level 6 in the second subsystem. 



6 



-1.4 



-1.8 



n=7-\^ 
n=6 \ 




(a) 





1.40114 



A 
A 



1.40117 



With increasing the level up to n = 15, we numer- 
ically make the RG analysis of the bicritical behavior. 
We first solve Eq. (27) and obtain the bicritical point 
(A* , i3* ) of level n and the pair of critical stability multi- 
pliers (At „, A2 p) of level n. Next, we use the formulas of 
Eqs. ( pll) and (^2|) and obtain the parameter and orbital 
scaling factors of level n, respectively. These numerical 
RG results for the first and second subsystems are listed 
in Tables IV and 0, respectively. Note that the accuracy 
in the numerical RG results is remarkably improved with 
the level n and their limit values agree well with those 
obtained by a direct numerical method. 



-1.174 



-1.184 



1 

n=6 ^"^^ 


' 1 


(b) 

1 





1.09009 



1.09010 



B 



FIG. 3. Plots of (a) the first stability multipliers Xi^n{A) 
versus A and (b) the second stability multipliers A2,,i(v4g,i3) 
versus B for the cases n = 6,7. In (a), the intersection point, 
denoted by the solid circle, of the two curves Ai,6 and Ai,7 
gives the point (Ag,Ai g) of level 6. As n 00, {A^,Xl „) 
converges to its limit point {A*,Xl). Similarly, in (b), the 
intersection point, denoted also by the solid circle, of the two 
successive curves A2, 6(^61 ^) ^^'^ A2, 7(^61 ^) gives the point 
(-Bg, A2.6) of level 6. As n — > 00, (Bi^, A2.„) also approaches 
its limit point (B*, A2). For other details, see the text. 



TABLE IV. Sequences of the critical point, the first critical stability multiplier, the parameter 



and orbital scaling factors, {^4^}, {A^ 


n}, {3i,n} and {ai,n}, in 


the first subsystem. For comparison. 


we 


also list the results obtained by a 


direct numerical method 


in the last row. 




n 


A*i 


Al,n 


Sl,n 




6 


1.401 155 189 088 9291 


-1.601 191211 1212 


4.669 203 072 1 


-2.502 620 459 5 


7 


1.401 155 189 092 133 2 


-1.601 191342 5171 


4.669 201428 5 


-2.502 845 988 3 


8 


1.401 155 189 092 048 4 


-1.601 191326 288 7 


4.669 2016314 


-2.502 894 652 


9 


1.401 155 189 092 050 7 


-1.601 191328 294 3 


4.669 201606 3 


-2.502 905 037 7 


10 


1.401 155 189 092 050 6 


-1.601 191328 046 4 


4.669 201609 4 


-2.502 907 267 8 


11 


1.401 155 189 092 050 6 


-1.601 191328 077 


4.669 201 609 1 


-2.502 907 744 9 


12 


1.401 155 189 092 050 6 


-1.601 191328 073 2 


4.669 201 609 1 


-2.502 907 847 2 


13 


1.401 155 189 092 050 6 


-1.601 191328 073 7 


4.669 201 609 1 


-2.502 907 8691 


14 


1.401 155 189 092 050 6 


-1.601 191328 073 6 


4.669 201 609 1 


-2.502 907 873 8 


15 


1.401 155 189 092 050 6 


-1.601 191328 073 6 


4.669 201 609 1 


-2.502 907 874 8 




1.401 155 189 092 050 6 


-1.601 191328 073 6 


4.669 2016091 


-2.502 907 875 1 
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TABLE V. Sequences of the critical point, the second crit- 
ical stability muhiplier, the parameter and orbital scaling fac- 
tors, {B*}, {A2,„}, {52,n} and {02, n}, in the second subsys- 
tem. For comparison, we also list the results obtained by a 
direct numerical method in the last row. 



n 


B^ 




52.n 


a2,n 


6 


1.090 092 490 313 


-1.177 467 


2.395 07 


-1.502 785 


7 


1.090 094 351702 


-1.178 671 


2.393 58 


-1.503173 


8 


1.090 094 321847 


-1.178 625 


2.393 59 


-1.504 426 


9 


1.090 094 328 376 


-1.178 649 


2.393 10 


-1.504 894 


10 


1.090 094 347 652 


-1.178 820 


2.392 80 


-1.504 993 


11 


1.090 094 348 817 


-1.178 844 


2.392 81 


-1.505 163 


12 


1.090 094 348 536 


-1.178 830 


2.392 78 


-1.505 263 


13 


1.090 094 348 675 


-1.178 847 


2.392 74 


-1.505 280 


14 


1.090 094 348 704 


-1.178 856 


2.392 73 


-1.505 296 


15 


1.090 094 348 701 


-1.178 853 


2.392 73 


-1.505 311 




1.090 094 348 701 


-1.178 85 


2.392 7 


-1.505 318 



IV. SUMMARY 

We have studied the scaUng behavior of period dou- 
blings near the bicritical point, corresponding to a thresh- 
old of chaos in both subsystems. For this bicritical case, a 
new type of non-Feigenbaum critical behavior appears in 
the second (response) subsystem, while the first (drive) 
subsystem is in the Feigenbaum critical state. Employ- 
ing the truncation and eigenvalue-matching methods, we 
made the RG analysis of the bicritical behavior. For 
the case of the truncation method, we analytically ob- 
tained the fixed point, associated with the bicritical be- 
havior, and its relevant eigenvalues. These analytic RG 
results are not bad as the lowest-order approximation. 
To improve the accuracy, we also employed the numer- 
ical eigenvalue-matching RG method, and obtained the 
bicritical point, the parameter and orbital scaling factors, 
and the critical stability multipliers. Note that the accu- 
racy in the numerical RG results is improved remarkably 
with increasing the level n. Consequently, these numer- 
ical RG results agree well with the results obtained by 
a direct numerical method. The results on the bicriti- 
cal behavior in the abstract system of unidirectionally- 
coupled ID maps are also confirmed in the real system 
of unidirectionally-coupled oscillators 11^1 . 
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